Dark Energy versus Modified Gravity 
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There is now strong observational evidence that the expansion of the universe is accelerating. The 
standard explanation invokes an unknown "dark energy" component. But such scenarios are faced 
with serious theoretical problems, which has led to increased interest in models where instead General 
Relativity is modified in a way that leads to the observed accelerated expansion. The question then 
arises whether the two scenarios can be distinguished. Here we show that this may not be so easy, 
demonstrating explicitely that a generalised dark energy model can match the growth rate of the 
DGP model and reproduce the 3+1 dimensional metric perturbations. Cosmological observations 
are then unable to distinguish the two cases. 
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INTRODUCTION 

The observed accelerated expansion of the late-time 
universe, as evidenced by a host of cosmological data like 
type la supernovae (SN-Ia) the cosmic microwave 
background radiation Q and large scale structure Q 
came as a great surprise to cosmologists. Although it 
is straightforward to explain the effect within the frame- 
work of Friedmann- Robertson- Walker cosmology by in- 
troducing a cosmological constant or a more general (dy- 
namical) dark energy component, all such explanations 
give rise to severe coincidence and fine-tuning problems. 

An alternative approach postulates that General Rel- 
ativity is only accurate on small scales and has to be 
modified on cosmological distances. This in turn leads to 
the observed late-time acceleration of the expansion of 
the universe [1, @, 0| ■ One of the best-studied exam- 
ples is the Dvali-Gabadadze-Porrati (DGP) brane-world 
model [8( , in which the gravity leaks off the 4-dimensional 
Minkowski brane into the 5-dimensional bulk Minkowski 
space-time. On small scales the gravity is bound to the 
4-dimensional brane and the Newtonian gravity is recov- 
ered to a good approximation. 

One important question is whether such a scenario can 
be distinguished from one invoking an invisible dark en- 
ergy component. It is well known that any expansion 
history (as parametrised by the Hubble parameter H{t)) 
can be generated by choosing a suitable equation of state 
for the dark energy (parameterised by the equation of 
state parameter w = p/p of the dark energy). This can 
for example be seen from Eq. (1) of Let us illustrate 
this explicitely for the DGP model, for which the Hubble 
parameter evolves as 
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where r c , the crossover scale, separates the 5D and the 
4D regimes. It has to be of the order of 1/Hq in or- 
der to generate late-time acceleration. Since matter is 



conserved on the brane, p m satisfies the usual conserva- 
tion equation. Comparing this to the normal Friedmann 
equation with an additional dark energy component, we 
see that we can move the crossover term to the right hand 
side and think of it as a dark energy contribution with 
Pde = 3H/(8irGr c ). Looking at the conservation equa- 
tion we find that it is solved if the effective dark energy 
has an equation of state given by 
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Consequently, it is impossible to rule out "dark energy" 
based on measurements of the cosmic expansion history 
(e.g. SN-Ia data). 

Recently there have been claims that it is instead pos- 
sible to use the growth rate of structures for this purpose 
[fil E3, EH [HI (but see also [IH for cautionary remarks). 
This is based on the observation that we can fix the equa- 
tion of state parameter uu of the dark energy from back- 
ground data and then predict the evolution of the dark 
matter perturbations in a standard cosmological model 
with dark energy. If the observed growth rate is different 
from the predictions, then general relativity with dark 
energy would be ruled out. 

However, in this paper we will show that this conclu- 
sions makes additional, very strong assumptions about 
the nature of the dark energy, and that in general the 
growth rate of structure is not sufficient to distinguish 
between dark energy models and modifications of grav- 
ity. We will show how the dark energy perturbations 
influence the dark matter and the metric perturbations, 
and provide an explicit example of a dark energy model 
which reproduces the 3+1 dimensional metric perturba- 
tions of the DGP scenario. 



SETTING THE STAGE 

We start by discussing the fluid perturbations in stan- 
dard 3+1 dimensional cosmologies. The perturbations in 
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the energy density are given by S = Sp/ p and to repre- 
sent the fluid velocity we use V = ikfT^j p. Working in 
the Newtonian (longitudinal) gauge, the metric can be 
written as 

ds 2 = - (1 + 24>) dt 2 + a 2 (1 - 2(f) dxidx 1 (3) 

with two scalar potentials <p and ip describing the per- 
turbations in the metric. Perturbations in cosmic fluids 
evolve according to [3] 
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where a prime denotes a derivative with respect to the 
scale factor a. The physical properties of the fluid are 
given by the anisotropic stress a and the pressure pertur- 
bation Sp (in general both can be functions of k). The 
latter is often parametrised in terms of the rest-frame 
sound speed c 2 , 



Sp = c 2 Sp + 3Ha(c 2 s - c 2 a )p— 



(6) 



where c 2 = p/p is the adiabatic sound speed. Collision- 
less cold dark matter has zero pressure (w rn = 0), van- 
ishing sound speed (c~ m = 0) and no anisotropic stress 
(a m — 0). For the dark energy all these quantities are a 
priori unknown functions and have to be measured. For 
the special case of dark energy due to a minimally cou- 
pled scalar field we have a variable w (corresponding to 
the choice of the scalar field potential, and fixed by the 
expansion history of the universe), c 2 DE = 1 and a = 
(see e.g. (3,0). 

The perturbations in different fluids are linked via the 
perturbations in the metric <f> and ip. Introducing the 
comoving density perturbation A = 6 + 3HaV/k 2 , their 
evolution in the standard cosmology is given by 

k 2 ^ = -4nGa 2l £2piAi (7) 

i 

k 2 ((j)-i{j) = 12nGa 2 ^2{l + w l )p l a l (8) 

i 

where the sum runs over matter and dark energy in our 
case. 

The quantity of interest to us is the growth factor 
g = A m /a which parameterises the growth of structure 
in the dark matter. The growth factor is normalised so 
that g = 1 for a <C 1 (using that A m oc a during matter 
domination and on sub- horizon scales). For dcfinitcncss 
we fix k = 200/f/o for the numerical results. We assume 
that g is an observable quantity (even though of course 
large scale structure surveys observe luminous baryonic 
matter, not dark matter, adding yet another layer of com- 
plications) . 
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FIG. 1: This figure shows how the growth of the matter 
perturbations depends on the clustering properties of the 
dark energy. From the top downward the sound speed is 
<? s — —2 x 10~ 4 (cyan dash-dotted line), (? s — — 10~ 4 (ma- 
genta long dashed line), c? s = (blue dotted line) and <? s = 1 
(red dashed line). For comparison we also plot the growth 
factor of the DGP model (black solid line) . 



THE IMPORTANCE OF DARK ENERGY 
PERTURBATIONS 



We start by noticing that the growth factor is not 
uniquely determined by the expansion history of the uni- 
verse (and hence wbe)- Although the main effect of the 
dark energy is to change H, leading to g < 1 at late times, 
there is an additional link through the gravitational po- 
tential ip. Different dark energy perturbations will lead 
to a different evolution of ip, which can modify the be- 
haviour of g. Conventionally one assumes that the dark 
energy perturbations are unimportant, e.g. 



Uj. This 



is a good assumption for scalar field dark energy where 
the high sound speed prevents clustering on basically all 
scales. However, a small sound speed Cj DE «0 is not ex- 
cluded. Indeed, it could even be negative, leading to very 
rapid growth of the dark energy perturbations. It could 
also vary in time. We show in Fig. Q] how the growth fac- 
tor of the dark matter changes in response to large dark 
energy perturbations [23| . 

What happens is that, as we decrease the sound speed, 
the dark energy is able to cluster more and more. The in- 
creased dark energy perturbations lead to enhanced met- 
ric perturbations. The dark matter in turn falls into the 
potential wells created by the dark energy, leading to an 
increase of the growth factor. Although clearly g is not 
uniquely determined by u>de, we notice that it always 
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increases as we decrease c 2 DE (at least as long as the lin- 
earised theory is applicable, see also [HI). Looking at the 
evolution equations ^ and (O for er = (<^> <j) = i/>) we 
see that the response of the fluids to the metric perturba- 
tions is governed by the sign of 1 + w. Non- phantom dark 
energy (as required to mimic the DGP expansion history) 
clusters therefore in fundamentally the same way as the 
dark matter and can only increase the growth of matter 
relative to the case of negligible dark energy perturba- 
tions (excluding highly fine-tuned initial conditions). 

So although the dark energy perturbations can influ- 
ence the growth factor of the dark matter, they only seem 
capable of enhancing it. But Fig.Q]also shows the predic- 
tion for the growth factor in the DGP model from [l9j |. 
and it is smaller than the one of a smooth dark energy 
component. We therefore need to change something else 
if we want to mimic DGP with dark energy. For this we 
need to take a closer look at the DGP model. 
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ANISOTROPIC STRESS AND MODIFIED 
GRAVITY MODELS 

An important aspect of DGP and other brane-world 
models is that the dark matter does not see the higher- 
dimensional aspects of the theory as it is bound to the 
three-dimensional brane. Its evolution is then the same 
as in the standard model. The modifications appear only 
in the gravitational sector, represented by the metric per- 
turbations. 

The metric perturbation in DGP can be written as 
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where the parameter (3 is defined as: 
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The dark matter does not care if the metric perturba- 
tions are generated (in addition to its own contribution) 
by a modification of gravity or by an additional dark en- 
ergy fluid. Its response to them is identical. Or to put it 
differently, if the dark energy and dark matter together 
can create the cf> and tp of Eqs. ((9]) and (fT0|) then the 
growth factor (and indeed all other cosmological observ- 
ables) will be the same as in the DGP scenario. 

We see immediately that in order to generate these 
metric perturbations we will need to introduce an 
anisotropic stress since 4> ^ tp. This seems to be a very 
generic property of modified gravity that is also present in 
f(R) models 2l| and has been noticed before. We plot in 



FIG. 2: In this figure we show how the anisotropic stress of 
the dark energy affects the growth of the dark matter pertur- 
bations. The red dashed line corresponds to scalar field dark 
energy with c 2 = 1 and a — 0. The dotted blue line shows 
how the dark matter growth factor decreases for a constant 
ode = —0.1. The long-dashed magenta line uses the theoreti- 
cal anisotropic stress of Eq. (112[) with c 2 s = 1, which suppresses 
the growth of the matter perturbations too much. Finally, the 
dash-dotted cyan line (nearly on top of black solid DGP line) 
uses the same ode but sets the pressure perturbation of the 
dark energy to Sp = (1 + w)pa in its rest frame. 



Fig. [2] again the growth factor for scalar field dark energy 
and the DGP model, but now also a family of dark en- 
ergy models with non-vanishing anisotropic stress a. We 
notice that these models can easily suppress the growth 
of perturbations in the dark matter for a < and mimic 
the behaviour of the DGP model. 

Formally we can recover the DGP metric perturbations 
by choosing 
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for the anisotropic stress of the dark energy, if we can 
also generate dark energy perturbations with 
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dp 
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We notice that these are very large dark energy pertur- 
bations. Indeed, if we keep c\ — 1 and set a to the 
expression (|L2[) we suppress the growth of the matter 
perturbations too much, see Fig. [2j Since (3 < the 
large dark energy perturbations of Eq. (fT3|) then increase 
the matter clustering back to the DGP value. 

The required size of the dark energy perturbations in 
itself is no problem, as we can lower the sound speed 
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and even make it negative. However, while for a — we 
were not able to decrease A m with the help of the dark 
energy perturbations, we find that with a large, nega- 
tive anisotropic stress we are unable to increase it. The 
required anisotropic stress is far larger than the gravi- 
tational potential ip, and it starts to be the dominant 
source of dark energy clustering in Eq. (JSJ). As it enters 
with the opposite sign it now leads to anti-clustering of 
the dark energy with respect to the dark matter which 
feels only ip (ie. dark matter overdensities are dark en- 
ergy voids). There is still enough freedom in the choice 
of a to match the growth factor very precisely, but if we 
could measure both <p and ip separately then we could 
detect the differences between the two models. 

Is it really not possible to match both ip and (p of 
the DGP model within a generalised fluid dark energy 
model? Yes, it is: The metric perturbations have two 
degrees of freedom, and we do have two degrees of free- 
dom of the dark energy to adjust, a and Sp. As it turns 
out, the parametrisation in terms of the rest-frame sound 
speed is too restrictive. This can happen for example if 
the dark energy is not composed of a single fluid, see 
e.g. the discussion in Allowing free use of the pres- 
sure perturbations, we can choose them for example to 
cancel the direct effect of a onto the dark energy pertur- 
bation in Eq. ([5]) by setting Sp = (l + w)pa. This reverses 
the sign of Ade, and minor adjustments to the pressure 
perturbations can then provide the required match to 
A m . For the cyan dash-dotted curve in Fig. [2] we set 
Sp = (1 + w)p<r + SHac^pV/k 2 , ie. we cancelled the con- 
tribution of the anisotropic stress in the dark energy rest 
frame. This provides a very good solution to Eqs. (fT^)) 
and (|13p during matter domination. It is easy to improve 
the solution to the point where it is impossible to distin- 
guish observationally between the DGP scenario and this 
generalised dark energy model. 

Is linear perturbation theory still valid with such a 
large anisotropic stress? Using Eq. Q13p we can rewrite 
Eq. (TT2]) as ctde = — 2/(3(1+wde))Ade- The anisotropic 
stress is therefore comparable in size to A m and Ade, and 
at high redshift DGP approaches GR. It is thus safe to 
study the dark energy with linear perturbation theory as 
long as the dark matter perturbations stay in the linear 
regime, even in the presence of the anisotropic stresses. 



CONCLUSIONS 

We have shown in this letter that the growth factor 
is not sufficient to distinguish between modified gravity 
and generalised dark energy, even if the expansion his- 
tory (and so the effective equation of state of the dark 
energy) has been fixed by observations. We have also 
demonstrated that in some cases (notably DGP) the dark 
energy can match the metric perturbations completely 
so that cosmological observations cannot distinguish be- 



tween the two possibilities. 

Although the construction of a matching dark energy 
model for the DGP case may seem very fine tuned, we 
are here more concerned with the question to what de- 
gree this is possible at all. Just measuring a growth factor 
that does not agree with scalar field dark energy is not 
sufficient to rule out "dark energy" and General Relativ- 
ity. But clearly, if the expansion history and the growth 
of matter perturbations were to match those predicted 
from a physically motivated and self-consistent modified 
gravity model, a statistical analysis would rule out a fine 
tuned dark energy model. However, we should not for- 
get that as observations seem to indicate wde w —1 it is 
rather the modified gravity models that are about to be 
ruled out [22| or look increasingly fine tuned. 

M.K. and D.S. are supported by the Swiss NSF. It is a 
pleasure to thank Chiara Caprini and Ruth Durrer for in- 
teresting discussions, and Eric Linder and Roy Maartens 
for comments on the draft. 



Martin.Kunz@physics.unige.ch 



[2 
[3 

K 
[5 

[6] 
[7] 

[8 

[9 

[10 

[11 
[12 

[13 

[14 

[is; 

[16 

[17 
[18 
[19 
[20 
[21 

[22 

[23 



Electronic address: 

Electronic address: Domenico.Sapone@physics.unige.ch 
A.G. Riess et al, Astron. J. 116, 1009 (1998); S. Perl- 
mutter et al, Astrophys. J. 517, 565 (1999). 

D. N. Spergel et al., ApJS 148, 175 (2003). 

M. Tegmark et al., Phys. Rev. D 69, 103501 (2004); M. 

Tegmark et al., Astrophys. J. 606, 702 (2004). 

C. Deffayet, Phys. Lett. B 502, 199 (2001) 

P. Binetruy, C. Deffayet, U. Ellewanger and D. Langlois; 

Phys. Lett. B 477 285 (2000). 

R. Maartens; |astro-ph/0602415| 

S. Capozziello, V.F. Cardone, S. Carloni and A. Troisi, 
Int. J. Mod. Phys. D 12, 1969 (2003). 
G.R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 
484, 112 (2000). 

C. Bonvin, R. Durrer and M. Kunz, Phys. Rev. Lett. 96, 
191302 (2006). 

E. V. Linder, Phys. Rev. D 72 043529 (2005). 

D. Huterer and E.V. Linder, |astro-p h/0608681 

D. Polarski, astro-ph/0605532 

E. V. Linder, astro-ph/0610173 

CP. Ma and E. Bertschinger, Astrophys. J. 455, 7 

(1995). 

W. Hu, |astro-ph/040 2060 

M. Kunz and D. Sapone, Phys. Rev. D 74, 123503 (2006). 

W.J. Percival, A&A 443 , 819 (2005). 

S. Dutta and I. Maor, |astn> ph/061202~7| 

K. Koyama, R. Maartens, JCAP 0601, 016 (2006). 

K. Koyama, JCAP 0603, 017 (2006). 

V.F. Mukhanov, H.A. Feldman and R.H. Brandenberger, 

Phys. Rep. 215, 206 (1992). 

R. Maartens and E. Majerotto, Phys. Rev. D 74, 023004 
(2006). 

We emphasize that we discuss only how the dark energy 
perturbations can modify the behaviour of the dark mat- 
ter, without taking into account limits from observations. 



